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Abstract 

We present a series of numerical experiments that model the evolution of magnetic 
flux tubes with a different amount of initial twist. As a result of calculations, tightly 
twisted tubes reveal a rapid two-step emergence to the atmosphere with a slight 
slowdown at the surface, while weakly twisted tubes show a slow two-step emergence 
waiting longer the secondary instability to be triggered. This picture of the two-step 
emergence is highly consistent with recent observations. These tubes show multiple 
magnetic domes above the surface, indicating that the secondary emergence is caused 
by interchange mode of magnetic buoyancy instability. As for the weakest twist case, 
the tube exhibits an elongated photospheric structure and never rises into the corona. 
The formation of the photospheric structure is due to inward magnetic tension force 
of the azimuthal field component of the rising flux tube (i.e., tube's twist). When 
the twist is weak, azimuthal field cannot hold the tube's coherency, and the tube 
extends laterally at the subadiabatic surface. In addition, we newly find that the 
total magnetic energy measured above the surface depends on the initial twist. Strong 
twist tubes follow the initial relation between the twist and the magnetic energy, while 
weak twist tubes deviates from this relation, because these tubes store their magnetic 
energy in the photospheric structures. 

Key words: magnetohydrodynamics: MHD, methods: numerical. Sun: corona. 
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1. Introduction 



Flux emergence is one of the key mechanisms in various solar activities. It is widely 
accepted that the emerging flux has a form of a twisted flux tube so as not to be collapsed 
by the convective motions during its ascent in the solar interior. Emerging flux transports 
magnetic energy and helicity from the convection zone to the atmosphere, which yields active 
regions including sunspots. Magnetic helicity in the corona stores free energy that can be 
released in the forms of flares and coronal mass ejections (CMEs) (e.g. Heyvaerts et al. 1977). 

Many numerical experiments have revealed the dynamics of the flux emergence. 
Schiissler (1979) carried out two-dimensional magnetohydrodynamic (MHD) simulations to 
study the cross- sectional evolution of the emerging flux tube (see also Moreno-Insertis & Emonet 
1996; Emonet & Moreno-Insertis 1998). Shibata et al. (1989) calculated the two-dimensional 
evolution of the undular mode of magnetic buoyancy instability (Parker instability: Parker 
1966) to reproduce the formation of an f2-shaped coronal loops. Toriumi & Yokoyama (2010) 
and Toriumi & Yokoyama (2011) gave numerical studies of emerging fluxes from much deeper 
convection zone (~ —20,000 km) to the corona. The three- dimensionality also exerts an influ- 
ence on emerging process of magnetic flux evolution. Matsumoto et al. (1993) produced the flrst 
three-dimensional work of the Parker instability using a magnetic flux sheet and a flux tube. 
Fan (2001) compared her numerical results of the twisted tube's emergence with observations 
of an active region. 

In this paper, we perform three-dimensional simulations of the twisted emerging flux 
tube from the uppermost convection zone to the corona. Our aim is to study the effect of the 
initial twist on the emergence process. A series of parametric studies on the flux tube's twist 
was done by Murray et al. (2006). Our work is dedicated to further detailed analyses, especially 
focusing on the effect of the initial twist on the resulting tube's structure (photospheric lateral 
expansion and multiple magnetic domes) and on the consequent coronal magnetic energy. 

For numerical experiments, we used the same conditions as those by Murray et al. (2006); 
we calculated ten cases of different twist parameters that cover their three runs. As a result of 
experiments, we found that the evolution depends on the initial twist. When the twist is strong 
enough, the evolution to the corona reveals two-step way, showing a deceleration and a lateral 
expansion near the solar surface, although the case with weaker twist spends more time waiting 
for the secondary emergence to occur (Magara 2001; Archontis et al. 2004; Murray et al. 2006). 
This picture of the two-step emergence is highly consistent with recent observations by Otsuji 
et al. (2010), especially its horizontally expanding speed and the rising speed. In addition to the 
confirmation of the results by Murray et al. (2006), it is also found that multiple magnetic domes 
are built and plasma accumulates in between the domes when the secondary emergence starts. 
At this moment, the direction of the field lines is almost perpendicular to the alignment of the 
domes, indicating that the second-step emergence is due to the interchange-mode instability. 
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If the initial twist is too weak, the tube extends widely near the surface and further evolution 
never takes place, because the magnetic tension force of the azimuthal component cannot hold 
the tube's coherency. 

Also, we newly found that the total magnetic energy measured above the surface relies 
on the initial twist. In the strong twist regime, the resulting magnetic energy follows the initial 
relation between the twist and the magnetic energy. In the weak twist regime, however, the 
magnetic energy deviates from the initial rule, because the tube with weak twist stores magnetic 
energy around the photosphere. 

The rest of the paper is organized as follows. In Section 2, we describe the numerical 
model. The simulation results are shown in Section 3. Summary and discussion are given in 
Section 4 and 5, respectively. 

2. Numerical Setup 

In numerical simulations, we solve nonlinear, time-dependent, compressible three- 
dimensional MHD equations. We take a rectangular computation box with three-dimensional 
Cartesian coordinates (x, z), where the 2-coordinate increases upward. The medium is as- 
sumed to be an inviscid perfect gas with a specific heat ratio 7 = 5/3. The basic equations in 
vector form follows: 
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where U is the internal energy per unit mass, I the unit tensor, the Boltzmann constant, 
m the mean molecular mass, and g = (0,0,— (70) the uniform gravitational acceleration. Other 
symbols have their usual meanings (p is for density, V velocity vector, p pressure, B magnetic 
field, c speed of light, E electric field, and T temperature). 
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To make above equations dimensionless, we introduce normalizing units of length Hq, 
velocity Cso, time tq = Hq/Cso, and density po, where Hq = k-QTo/ (mgo) is the pressure scale 
height, Cso the sound speed, and po the density at the photosphere, respectively. The gas 
pressure, temperature, magnetic field strength, and energy are normalized by the combinations 
of the units above, i.e., po = PqC^q, Tq = mC^Q/{'ykB), Bq = (poC^^o)^^^' and Eq = pqC^qH^, 
respectively. The gravity is given as go = C'^q/ (7-^0) by definition. For comparison of numerical 
results with observations, we use i^o = 170 km, Cso = 6.8 km s~^, tq = Hq/Csq = 25 s, and 
Po = 1-4 X 10"^ g cm""^, which are typical values for the solar photosphere. Then, po = 6.3 x 
10^ dyn cm^^ Tq = 5600 K, Bq = 250 G, and Eq = 3.1 x 10^6 erg. 

The initial background stratification consists of three regions: an adiabatically stratified 
convective layer, a cool isothermal photosphere/chromosphere (afterward, we simply call it 
photosphere), and a hot isothermal corona. The photosphere and the corona are smoothly 
connected by the transition region. We take z / Hq = to he the base height of the photosphere, 
and the starting height of the transition region and the corona are Ztr/ Ho = 10 and Zcor/ Hq = 20, 
respectively. The initial temperature of the photosphere and the corona are Tph/To = 1 and 
Tcot/Tq = 150, respectively. The background temperature distribution in the transition region 
is given by 

{z — Ztv)/{ZcoT—Ztr) 

(8) 
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is the adiabatic temperature gradient (subscript s is for surrounding distribution). The initial 
gas pressure and density profiles are defined by solving one- dimensional hydrostatic equation 

-^p,{z)+p,{z)go = (11) 

and the equation of state (7) on the basis of the temperature distribution above. 

The initial magnetic fiux tube is embedded in the convection zone at Ztnhc/Ho = —10. 
The longitudinal and azimuthal component of the fiux tube are described as follows: for a 
radial distance from the axis r = [{y — ytubc)^ + (-2 ~ ^tubc)^]"'^''^. 
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where (ytubc^tube) = (0, — lOi^o) is the tube center, -Rtubc the radius, q the twist parameter, and 
Btube the magnetic field strength at the axis. We take Rtuhc/Ho = 2.5 and Btuhc/Bo = 15, i.e., 
these parameters are almost the same as those of Murray et al. (2006). For pressure balance 
between the fiux tube and the surrounding medium, the gas pressure inside the tube is obtained 
as Pi = Ps + 5pcxc, where 



and X/Hq = 20. That is, the fiux tube is most buoyant at the middle of the tube {x / Hq = 0) , 
and the buoyancy diminishes as |x| increases. 

Here, we investigate ten parameters of q, which are qHo = 0.5, 0.4, 0.3, 0.25, 0.2, 0.175, 
0.15, 0.125, 0.1, and 0.05. The plasma beta (/3 = Sttp/B"^) at the tube center is /3 ~ 3 at the 
initial state. The initial background stratification (gas pressure, density, and temperature) and 
the magnetic pressure along x/Hq = h/Hq = of the case qHo = 0.2 are indicated in Figure 1. 

The simulation domain is taken as (—120, —120, —20) < {x/Hq, y/Ho, z/Hq) < 
(120, 120, 150), resolved by 256 x 256 x 256 grids. The grid spacings for x, y, and z directions 
are Ax/Hq = Ay/Ho = 0.5 for (-40, -40) < {x/Hq, y/Ho) < (40, 40), and Az/Hq = 0.2 for 
—20 < z/Hq < 20, respectively. Outside this range, the mesh sizes gradually increase. We as- 
sume periodic boundaries for horizontal directions and symmetric for vertical. A wave-damping 
region is attached near the top boundary. We use the modified Lax-Wendroff scheme version 
of the CANS (Coordinated Astronomical Numerical Software) code (see Toriumi & Yokoyama 
2010). 

3. Results 

3.1. Overview of the Results 

Figure 2 shows the time evolution of the fiux tube with the twist qHo = 0.2. In this 
Figure, we plot the logarithmic field strength log (|i?|/-Bo) and the photospheric magnetogram 
Bz/Bq. In each panel, the region x/Hq < and y/Ho>0 is shown. Initially, the fiux tube is 
embedded at ^tubcZ-f^o = ~10, and is slightly buoyant around the tube center —20 < x/Hq < 20. 
The tube rises through the convection zone by magnetic buoyancy, and reaches the surface 
at t/ro = 20 (Figure 2(b)), while the outskirt of the rising portion {x/Hq ~ —20) begins to 
sink, because the fiuid is drained along the field lines from the apex of the rising tube. Due 
to the isothermal (i.e., strongly-subadiabatic) photosphere, the tube is decelerated and ex- 
pands laterally near the surface to make a "photospheric tongue" (Figure 2(c)). That is, the 
convectively stable photosphere inhibits an upward motion of the fiuid, and thus, the mag- 
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The density inside the tube is also defined as Pi = Ps + Spc^c, where 
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netic field cannot penetrate the photosphere only to escape in the horizontal direction. (It 
should be noted that the term "tongue" here is different from that used in Li et al. (2007) and 
Archontis & Hood (2010). The horizontal extension of the magnetic field is also observed in 
a recent radiative MHD calculation by Cheung et al. (2010). The lateral expansion speed is 
\Vy\ ~ 0.4Cso = 2.7 km s~^, i.e., a fraction of the photospheric sound speed. At this moment, 
the photospheric magnetogram shows a north-south (y-directional) magnetic distribution, and 
the total field strength is \B\ ~ 2Bo = 500 G and plasma beta is (3 = p/pmag ~ 2 around the 
photosphere. 

As the magnetic pressure gradient enhances, the second-step emergence takes place. 
In Figure 2(d), multiple expansions are observed: the rise velocity is about (0.3 — 0.5)Cso = 
2.0 — 3.4 km s~^. This multi-dome structure is more noticeable in a weaker twist case. Figure 
3(a) shows the magnetic field structure of the case qHo = 0.15 at t/ro = 80. The corresponding 
density structure and velocity vectors in the y/HQ = plane, and the field lines are indicated 
in Figure 3(b). From these figures, one can see four magnetic domes are built and the fluid 
is accumulated between the domes. In this region, the field lines are generally directed in the 
^/-direction. Therefore, this situation can be explained as a consequence of the interchange 
mode of the magnetic buoyancy instability, i.e., the wavenumber vector is perpendicular to the 
field lines. 

As time goes on, the flux tube expands both vertically and horizontally, while the pho- 
tospheric tongue also continues to expand laterally (Figure 2(e)). Finally, the flux tube makes 
a single dome of the height z/ifo ~ 60 on the pancake-like structure at the surface (Figure 2(f)). 
At the same time, the sunk part approaches the bottom of the simulation domain [z/ H^r^ —20). 

Overall evolution described above is similar to the observation by Otsuji et al. (2010). 
They found the lateral expansion with the speed of 2.9 km s~^ at the surface before further 
evolution occurred, which is consistent with our results of \Vy \ = 2.7 km s~^. The gradual rise 
speed of the secondary emergence was observed to be 2.1 km s~^, which is also consistent with 
our results of 2.0 — 3.4 km s~^. 

3.2. Parameter Study on the Twist Strength 

Figure 4 shows the height-time relation of the top of the tube for various twist cases. 
Here, in this figure, we plot the height at the highest portion of the emerging flux tube 
(^apex/-f^o)- The evolutions are found to depend on the initial twist, and each line distributes in 
a continuous fashion. It can be seen from this figure that almost all the tubes (0.5 > qHo > 0.1) 
show the two-step emergence to the corona. Although the rise times within the convection zone 
are similar to each other, tubes with weaker twists spend more time in the surface waiting for 
the second-step emergence to be triggered. The tube with qHo = 0.1 shows only a slight emer- 
gence in the atmosphere {z/Hq < 20). As for the weakest twist case with qHo = 0.05, further 
evolution never takes place within the elapse calculated (failed emergence). These results are 
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consistent with those of Magara (2001), Murray et aL (2006), and Toriumi & Yokoyama (2011). 

Cross-sections at x/Hq = plane of eight out of ten flux tubes are shown in Figure 5. 
These tubes are those who reach z/Hq = 40, and each figure shows the arrival at that height. 
As the initial twist qIIq becomes smaller, the lateral expansion at the surface (tongue-like 
structure around z/i^o ~ 0) is increasingly remarkable, because the more intense initial twist 
of the magnetic flux tube yields the stronger azimuthal magnetic tension force, and thus keeps 
the tube coherent. Therefore, the coronal field intensity also reduces with decreasing qHo, and 
the coronal structure with a weaker twist is more fragmented compared to the stronger twist 
cases. 

3. 3. Comparison of the Magnetic Field Structure at the Surface 

In this subsection, we compare the magnetic field structures of the different emergence 
cases to study the horizontal expansion and the mechanism of the second-step evolution at 
the solar surface. Figure 6 shows (top) the cross- sectional configuration of the flux tubes at 
x/Ho = with velocity vectors, {middle) the horizontal components of forces along the horizontal 
axis x/Hq = z/Hq = 0, and {bottom) the vertical components of forces along the vertical axis 
x/Hq = y / Hq = for cases with qEq = 0.4 (rapid emergence), 0.15 (slow emergence), and 0.05 
(failed emergence) at the time t/ro = 40. 

As can be seen from the top figures, the second-step evolution has already begun at this 
time for the case with qHq = 0.4, while the vertical expansion cannot be seen in qHq = 0.05 
case. It should be noted that the downfiow in the uppermost areas in Figures 6(b) and (c) 
is a reflected wave from the top boundary. Horizontally, for a strongest twist case (Figure 
6(d)), an inward magnetic tension force is dominant, and thus the inward total force keeps the 
tube from a lateral fragmentation. As the twist decreases, the magnetic tension reduces so 
that the total force is outward in a wider range for qHo = 0.05 (Figure 6(f)), resulting in the 
tube's fragmentation and the further expansion never to occur. As for vertical forces, magnetic 
pressure gradient is principal at front of the tube (Figure 6(g): 3 < z/Hq < 9). However, 
the total force is about zero within this area. That is, the second-step expansion is caused 
by the tube's magnetic pressure, while the tube is almost in a hydrostatic equilibrium with 
surrounding materials. For a weak twist case (Figure 6(i)), the magnetic pressure gradient is 
much less effective. Therefore, further rise cannot occur. 

3.4- Undulating Configuration of the Photospheric Field Lines 

Figure 7 shows the surface magnetogram Bz/Bq and the field lines above the surface for 
the medium twist case {qHq = 0.15) at the time t/ro = 125, namely, in the later phase. We 
confirm that some undulating field lines connect magnetic patches at the surface, and that, as 
time goes on, they gradually rise into the corona by forming longer fields. At this time, both 
near-surface undulating fields and coronal fields are directed almost parallel to the axis of the 
original flux tube, which is in contrast to the perpendicular fields observed in the earlier phase 
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of the emergence (see Figure 3(b)). This result is reminiscent of the "sea-serpent" field lines 
and the resistive emergence by Pariat et al. (2004). Undulating photospheric fields of a weakly 
twisted tube are also found by Archontis & Hood (2010). We will discuss this again in Section 



3.5. Magnetic Energy in the Atmosphere and the Initial Tube's Twist 

Figure 8 shows the initial twist qH^ and the total magnetic energy measured above the 
solar surface E^^^/ Eq when each tube arrives at z/Hq = AQ (see Figure 5). The magnetic energy 
above the surface is defined as: 



When the initial twist is large (gi^o > 0.2), the magnetic energy -Emag is found to obey a 
law. From the field configuration (12) and (13), the initial magnetic energy per unit volume 
can be calculated as 



that is, the initial tube's magnetic energy depends on q^. When the twist is strong, the az- 
imuthal magnetic tension force is more effective and the lateral expansion around the photo- 
sphere is less efficient that the initial magnetic energy is directly transported into the atmo- 
sphere above the surface. As a result of this regime, the observed magnetic energy within the 
atmosphere relies on a q^ function. 

Contrary to this, when the initial twist is weak {qH^ < 0.2), the consequent magnetic 
energy deviates from the q^ line and is negatively-correlated with the initial twist. It is specu- 
lated that, for these tubes with weaker twits, it takes a longer time to reach z/ Hq = 40 from 
the surface. It is because the photospheric field needs more time to satisfy the condition for 
the second-step emergence (i.e. Acheson 1979), for the weak azimuthal field causes the tube 
to expand horizontally around the surface. Since the magnetic energy is continuously trans- 
ported from below and the horizontal expanding velocity is almost the same for qH^ < 0.2 cases 
(Vy/Cso ~ 0.1), the magnetic energy of the photospheric tongue is expected to depend on the 
time lag between the tube's arrival at the photosphere and at z/Hq = 40. 

It can be concluded that the magnetic energy in the atmosphere correlates critically with 
the initial twist. In the strong twist regime, the energy-twist relation follows the initial q^ rule, 
because the tightly twisted tube does not exhibit a significant expansion near the photosphere. 
As for the weak twist cases, they depend on the time lag between reaching the photosphere and 
the corona. It is because the weaker twist tube takes more time to rise further, and, therefore, 
more magnetic energy is stored in the photospheric tongue. 



5. 




(16) 




(17) 
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4. Summary 



In this paper, we carried out three-dimensional MHD simulations to investigate the effect 
of the initial twist on the flux tube evolution. Here, we summarize the results: 

• Initially, the flux tube rises through the convection zone due to its magnetic buoyancy. 
Reaching the surface, the tube expands laterally to make a "photospheric tongue." The 
secondary emergence occurs after sufficient flux accumulates within the photosphere. Due 
to the interchange mode instability, the tube builds multiple domes above the surface, 
between which the fluids piles up. Finally, the flux tube arrives at z/iJo ~ 60 as a single 
dome. The overall emergence is consistent with the recent observations (e.g. Otsuji et al. 
2010). 

• We run ten twist cases to investigate the effect of the initial twist. Nine out of ten reach 
the coronal height {z/Hq > 20) showing two-step emergence, while the weakest twist case 
fails to rise further above the surface {qHo = 0.05). In the two-step emergence regime, 
the rise time becomes shorter with increasing initial twist, which is consistent with the 
previous calculations by Murray et al. (2006). The photospheric tongue is more noticeable 
in weaker twist case. 

• We study the force components at the solar surface for different twist cases at the time 
I/tq = 40. The stronger the initial twist is, the larger the inward magnetic tension is, 
resulting the tube keeps its coherency. In the weak twist case, the magnetic tension is 
much less effective, causing the tube distorted. At the same time, the strong twist tube 
rises further into the atmosphere mainly by the magnetic pressure gradient. 

• We found that the photospheric fields of the middle twist case {qHo = 0.15) undulate in 
the later phase of the emergence. The field lines gradually rise into the corona as longer 
loops. The photospheric and coronal fields are almost parallel to the axis of the initial flux 
tube. These features remind us of the resistive emergence model by Pariat et al. (2004). 

• We measure the magnetic energy -Emag above the surface. The energy plot follows the initial 

law when the twist is strong {qHo > 0.2), while, for weaker twist cases {qHo < 0.2), the 
energy depends on the time difference between reaching the surface and the corona. That 
is, weakly twisted tube takes more time for magnetic flux to accumulate near the surface 
and the secondary instability to be triggered. 

5. Discussion 

In Section 3, we showed the time-evolution of the twisted flux tube. When the second- 
step emergence starts, multiple domes are observed above the surface and fluid is trapped 
between the expanding magnetic structures (see Figure 3). At this time, field lines are directed 
perpendicular to the alignment of the magnetic domes. For middle twist tube, we also found 
undulating fields near the surface, emerging into the corona (see Figure 7). In this section, we 
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discuss these features in connection with future observations. 

5.1. Twist Intensity and the Interchange Instability 

In Section 3.1, we saw that, as the twist decreases, the interchange mode instabihty 
becomes more noticeable. However, this is contrary to the expectation that the azimuthal field 
should be less pronounced in a weaker twist case. It may be because, in a weak twist case, the 
tube extends laterally near the photosphere and thus the twist increases. As the tube develops 
the interchange instability, field lines perpendicular to the alignment of the magnetic domes 
become more pronounced (see Sections 3.2 and 3.3). 

5.2. Twist of the Actual Flux Tube in the Sun 

Multiple magnetic structures and the density accumulation between them are also found 
in previous observations and calculations. Pariat et al. (2004) found that photospheric fields are 
undulating at its earlier phase of the flux emergence event, and proposed a resistive emergence 
model that undulating multiple loops reconnect with each other to make larger coronal fields. 
Isobe et al. (2007) carried out two-dimensional MHD simulation to study the evolution of 
the serpentine magnetic loops (resistive emergence model), finding that density accumulates 
in between the magnetic loops; their elongated vertical plasma structures are similar to our 
results (for three-dimensional study, see Archontis & Hood (2009)). However, in our model, 
field lines are directed almost perpendicular to the alignment of the domes, which is against the 
observations of the undular field lines (e.g. Pariat et al. 2004). The direction of the field lines 
are the consequence of the initial tube's twist. Therefore, the difference between the present 
calculations and the observations indicates that the actual twist of the fiux tube beneath the 
surface may be much weaker than those assumed in our models (e.g. qHo = 0.15 for Figure 3). 

On the other hand, fiux tube with insufficient twist was found to fail to rise through the 
convection zone (Moreno- Insert is & Emonet 1996; Emonet & Moreno-Insertis 1998; Toriumi 
& Yokoyama 2011). It is because the weak azimuthal field of the fiux tube cannot hold its 
coherency during its ascent within the solar interior. Therefore, one of the important problems 
to be solved is the emergence of the fiux tube with much weaker twist {qHo < 0.1). 

We also found that the medium twist tube {qHo = 0.15) reveals the undulating fields at 
the surface, which gradually rise into the upper atmosphere as longer coronal loops (Section 7). 
These fields are directed parallel to the main axis of the initial flux tube. This picture seems 
well accorded with the resistive emergence model. However, it is in the later phase that this 
undulatory evolution is observed, and, in the earlier phase, the field lines are perpendicular 
to the original tube's axis (see Figure 3). Therefore, the reproduction of the undulating fields 
parallel to the original axis (Pariat et al. 2004) is not achieved. 
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5.3. Formation of Undulating Photospheric Fields 

Recently, Cheung et al. (2010) have conducted a radiative MHD simulation of the forma- 
tion of an active region. They showed that the rising flux tube flattens to make a pancake-like 
structure near the surface, and that the convective flows create serpentine field lines. 

In the present study, which does not include the convection, we also found the sideway 
expansion of the rising field at the photosphere (tongue) and the undulation of the photospheric 
fields in the later phase of the emergence of the weaker twist tube. These features are also 
confirmed by Archontis & Hood (2010); their simulations do not take account of the convection 
effects, either. Therefore, we can see that other mechanisms, apart from the convection, could 
also explain the formation of serpentine fields at the photosphere. 

5.4- For Future Observations 

In this paper, we found some aspects of the flux emergence event. One is the pho- 
tospheric tongue, i.e., the magnetic structure extending horizontally around the surface just 
before further evolution takes place. Temporally- and spacially-resolved spectroscopic observa- 
tions of the earlier phase of the flux emergence are required to study this magnetic extension 
at the photosphere. 

At the same time, we found that the initial twist of qHq ~ 0.1 at —1700 km is too strong 
to match the observations (as mentioned above in Section 5.2). Local and global helioseismology 
are needed to reveal the flux emergence (especially on the twist evolution) within the convection 
zone. The key issue is how weakly twisted flux tubes manage to rise through the solar interior. 

Numerical computations were carried out on NEC SX-9 at the Center for Computational 
Astrophysics, CfCA, of the National Astronomical Observatory of Japan, and on M System 
(Fujitsu FXl) of JAXA Supercomputer System. The page charge of this paper is partly sup- 
ported by CfCA. S. T. and T. Y. thank Dr. Y. Fan of the High Altitude Observatory, the 
National Center for Atmospheric Research. We thank the referee for helpful suggestions for 
improvements of this paper. 
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Fig. 1. One-dimensional (z-)distributions of the initial background density (thick solid line), pressure 
(dotted line), and temperature (dashed Hne). The magnetic pressure Pm = B'^ / (Stt) of the tube qH^ = 0.2 
along the vertical axis x/Hq = y/Ho = is overplotted with a thin solid line. 
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Fig. 2. Time evolution of tlie flux tube with the twist qHo = 0.2. Logarithmic field strength Iog(|i3|/i3o) 
and photospheric magnetogram B^/Bq are plotted. In each panel, the region x/Hq < and y/Ho > is 
shown. This figure is also available as an avi animation in the electronic edition. 
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Fig. 3. The flux tube with qH^ = 0.15 at the time t/ro = 80. (a) Magnetic field structure Iog(|i3|/i3o) and 
photospheric magnetogram Bz/Bq. Plotted colors are the same as those of Figure 2. (b) The corresponding 
logarithmic density profile \og{p/po) at y/Ho = plane (—40 < x/Hq < —5 and 5 < z/Hq < 17.5) with 
velocity vectors (white arrows), magnetogram at z/Hq = 5, and field lines (blue lines) are shown. 
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Fig. 4. Height-time relations for various twist cases. From left to right, each line shows the evolution of 
qHo = 0.5, 0.4, 0.3, 0.25, 0.2, 0.175, 0.15, 0.125, 0.1, and 0.05. It takes more time to rise into the corona 
as the twist strength decreases. 
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Fig. 5. Cross-sections at x/Hq = plane of eiglit out of ten flux tubes that reach z/Hq = 40. Each figure 
shows the logarithmic field strength log(|i?|/i?o) when the tube reaches that height. 
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Fig. 6. (Top) Cross-sectional configuration of the flux tubes at x/i?o = with velocity vectors, (middle) 
horizontal components of each force along the horizontal axis x/Hq = z/Hq = 0, and (bottom) vertical 
components of forces along the vertical axis x/Hq = y/Ho = for cases with qHo = 0.4 (rapid emergence), 
0.15 (slow), and 0.05 (failed) at the time t/ro ~ 40. Plotted lines are the total force (thick solid line), 
gas pressure gradient (dotted), magnetic pressure gradient (dashed), magnetic tension (dash-dotted), and 
gravity (thin solid), respectively. 
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Fig. 7. Surface magnetogram B^/Bq (color contour) and the field lines (blue lines) for the middle twist 
case {qHo = 0.15) at the time t/ro = 125 are shown. Note that some filed lines undulate near the surface 
and gradually rise into the corona as longer loops. 
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Fig. 8. Total magnetic energy above the surface -Emag = Jz>o-^^ /(^^) versus the initial twist q, 
measured at the time when each tube reaches z/Ho = 40. Quadratic and linear lines are overplotted with 
sohd lines. 
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